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Abstract. Let ^ H (x,t) = 2H /* 5(Bf - x)s 2H ~ 1 ds be the weighted local time of 
fractional Brownian motion B H with Hurst index 1/2 < H < 1. In this paper, we use 
Young integration to study the integral of determinate functions f R f(x)J? H (dx,t). 
As an application, we investigate the weighted quadratic covariation [f(B H ),B H ] < - w * > 
defined by 

n-l 

[f(B H ),B H ]r := Em 2Jf ^ ) - - B? k ), 

k=0 

where the limit is uniform in probability and t& = kt/n. We show that it exists and 

[/(B H ),BT' = - / f(x)J? H (dx,t), 
Jm 

provided / is of bounded p- variation with 1 < p < ^_ H - Moreover, we extend this 
result to the time-dependent case. These allow us to write the fractional Ito formula 
for new classes of functions. 

1. Introduction 

The fractional Ito's formula has played a central role in stochastic analysis for frac- 
tional Brownian motion and almost all aspects of its applications. But, the restriction 
of Ito's formula to functions with twice differentiability often encounter difficulties in 
applications. Extensions to less smooth functions are useful in studying many prob- 
lems such as partial differential equations with some singularities and mathematics of 
finance. For the classical Ito's formula, many authors have written extended versions of 
the Ito formula in order to relax this regularity condition. Still there are always some 
new problems requiring the use of "Ito formula" under lighter conditions. One approach 
of extending Ito's formula is by using local time-space calculus to absolutely contin- 
uous function F with locally bounded measurable derivative F' (see Bouleau-Yor [3], 
Eisenbaum [HE], Feng-Zhao [7], Follmer et al jSj, Peskir [23], Russo-Vallois [M], Yan- 
Yang [26], and the references therein). Moreover, the backward integral and quadratic 
covariation are fundamental tools in these discussions. However, all these extensions 
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are not effective on fractional Brownian motion B H since the quadratic covariation 
[f(B H ),B H ] of f(B H ) and B H satisfies 

[+00, it < H < ± 

for all t > 0, in general, where x 1— > /(cc) is a determinate function. As well-known, the 
fractional Brownian motion (fBm) with Hurst index H £ (0, 1) is a mean zero Gaussian 
process B H = {Bp, < t < T} with the covariance function 

E [B t H B?] = \ [t 2H + s 2H -\t- s\ 2H ] 

for all t, s ^ 0. For H = 1/2, B H coincides with the standard Brownian motion B. 
B H is neither a semimartingale nor a Markov process unless H = 1/2, so many of 
the powerful techniques from stochastic analysis are not available when dealing with 
B . However, as a Gaussian process, one have constructed the stochastic calculus of 
variations with respect to B H (see Biagini et al [2] and references therein). 

In order to motivate our subject, let us first recall some known results concerning 
the quadratic variation and Ito's formula. Let F be an absolutely continuous function 
with locally square integrable derivative /, that is, 



F(x) = F(0)+ f f(y)dy 
Jo 





with / being locally square integrable. Follmer et al [8] introduced the following Ito's 
formula: 

(1.1) F(B t ) = F(0) + f(B s )dB s + \ [f(B), B] t . 

If / £ C^IR), (jl.ip is the classical ltd formula. Moreover, the result has been extended 
to some semimartingales and smooth nondegenerate martingales (see Russo-Vallois [23] 
and Moret-Nualart [20]). However, it is obvious that the formula (jl.ip is not true for 
fBm B H if the corresponding stochastic integral J Q f(B^)dB^ is of Wick-Ito type 
and / C 1 (R). Thus, it is natural to ask whether the similar ltd formula for fBm 
B holds or not. This motivates the subject matters of this paper. Here, we only 
consider the case ^ < H < 1. This means that fBm B H is of finite quadratic variation 
provided | < H < 1. Recall that a process X is said to be of finite quadratic variation 
if quadratic variation [X, X] is finite. For any continuous finite quadratic variation 
process X we have (see, for example, Gradinaru et al [TT] ) 

(1.2) F(X t ) = F(0) + |* f(X s )d-X s + \ [f(X),X] t 

provided F G C 2 (]R) and F' = /, where the integral J* f{X s )d~ X s is the forward 
(pathwise) integral defined by 

f f(X s )d-X s = lim - f f{X s ){X s+£ - X s )ds 
Jo £ i° £ Jo 

and [f(X),X] t = f* f'(X s )d[X, X] s , and more precisely, for fBm B H with \ < H < 1 
we have 

f f(B»)d-B»= ff{Bf)dB? + H f f{B?)s 2H - l ds, 
Jo Jo Jo 
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where the integral J Q f(B^)dB^ is of Wick-Ito type. However, the formula (| 1 . 2 j) is 
only effective on twice-differentiable functions. It is impossible to list here all the 
contributors in previous topics. Some surveys and complete literatures could be found 
in Nualart [21], Biagini et al [2j, Hu [14J, Mishura [18], Russo-Vallois [25] and Gradinaru 
et al pa QT]. 

In the present paper, our aim is devoted to find a substitution tool of the quadratic 
variation such that Ito's formula similar to (II. lh holds for fBm with ^ < H < 1 whatever 
/ S C 1 (]R). For fBm, in [11] (see also [12] and the references therein) Gradinaru et 
al have introduced some substitution tools and studied some fine problems. They 
introduced firstly an ltd formula with respect to a symmetric-Stratonovich integral, 
which is closer to the spirit of Riemann sums limits, and defined a class of high order 
integrals having an interest by themselves. They also treated an Ito formula with 
respect to somehow any symmetric integral, introducing a large class of symmetric 
integrals via regularization, and fractional Brownian motion is not the only process for 
which their Ito formula is valid; there are easy extensions to a more general class of 
processes. However, our substitution of quadratic covariation and technique used here 
are different from theirs. Moreover, we use the integral of determinate functions with 
respect to the weighted local time of fBm B H 

[ [ g(x,s)^ H (dx,ds) 
Jo Jr 

and 

/ f(x)J? H (dx,t) 
Jr 

to weaken the hypothesis in some formulas, where ^ H (x,t) = 2H £ 5{Bf -x)s 2H ~ 1 ds 
is the weighted local time of fBm B . Though our method is only effective on fBm 
with I < H < 1, the merit here has been to concentration fully on fBm in order to get 
a stronger statement by fully using fBm's regularity. 

Recently, in [10] , [16] and [T7] Nourdin et al studied weighted power variations of fBm 
and introduced some fine results. More precisely, they showed that the convergence 

n 2H ^ 1 ^2h(B^ /n )(Bf k+1)/n -B^ /n ) — j^/ h{B»)ds a.s. 
k=0 J ° 

holds for all < H < 1, as n tends to infinity, where p > 1 is an integer number, h E 
C(R) and \xi v is the 2p-moment of a standard Gaussian random variable G ~ N(0, 1). 
As a corollary, the following convergence holds: 

(1.3) £ fc 2 ^~ W/J K + i)/n " Kn) — / HB^s^ds a.s. 
k=0 Jo 

for all < H < 1. This inspirits us to consider the weighted quadratic covariation 
[f(B H ),B H ]( w ) of f{B H ) and B H , defined by 

n-1 
k=0 
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where the limit is uniform in probability and = kt/n. Clearly, for / £ C (R) the 
convergence (|1.3|) implies that [f(B H ), B H ][ W ^ exists and equals to 



2H f f'(Bf)s^ds, 
Jo 
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and moreover the fractional Ito formula can be rewritten as 

F(2?f ) = F(0) + f f{B?)dB* + H f f(B?)s 2H - l ds 
Jo Jo 

^ F(0) + jf f(B?)dB? + i [/(-B^), B H ] ^ 

for all i G [0, T] and all < H < 1. This is the start point in this paper and we 
only consider fBm with i < H < 1. Our aims are to study the integrals with respect 
to the weighted local time of fBm B H and the weighted quadratic covariation, and to 
extend the above the fractional Ito formula to / G" C 1 (M). The case < H < \ will be 
discussed in a forthcoming paper. 

This paper is organized as follows. In Section [2] we present some preliminaries for 
fBm. To use Young integration to establish integral with respect to local time, in 
Section [3l we first investigate the power variation of the function x \— > J§f (x, t) for 
every t > 0. We show that 

lim J2\Sf H (ai +1 ,t)-Sf H (ai,t)\^ 



A 



exists in L 1 for every < t < T and \ < H < 1, where (A n ) is a partition of 
the interval [a, b] such that |A n | — > as n — > oo. Thus, as a direct consequence we 
construct a one parameter integral with respect to these local times with ^ < H < 1. 
The weighted quadratic covariation is considered in Section [H We show that if the 
measurable function x h- > /(x) is of bounded p- variation with 1 < p < and 

| < il < 1, then the weighted quadratic covariation B^l exists in L 1 and 

[/( J B // ), J B i/ ]^ ) = - / f(x)J? H (dx,t), t E [0,T], 
Jk 

and moreover the Ito formula 

F(Bf) = F(0) + J*f(B*)dB* + \ [f(B«),B»] { p 

holds for all absolutely continuous function F(x) = F(0) + Jq f(y)dy. In Section [5] we 
extend these results to the time-dependent case. 



2. Preliminaries on fractional Brownian motion 

In this section, we briefly recall some basic definitions and results of fBm. For more 
aspects on these materials we refer to Nualart [21J, Hu |14j . Biagini et al [2], Mishura [18] 
and the references therein. Throughout this paper we assume that \ < H < 1 is 
arbitrary but fixed and let B H = {Bf, < t < T} be a one-dimensional fBm with Hurst 
index H defined on P). Let (S)* be the Hida space of stochastic distributions 
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and let o denote the Wick product on (S)*. Then t i— ► Bp is differentiable in (5)*. 
Denote 

w m = d_B[_ g {sy _ 

We call the fractional white noise. For u : R + — > (5)*, in a white noise setting 

we define its (generalized) fractional stochastic integral of ltd type with respect to B H 
by 

(2.1) / u s dBf := f u s oW^ H) ds, 



whenever the last integral exists as an integral in (5)*. We call these fractional Ito 
integrals, because these integrals share many properties of the classical Ito integral. 
For any F £ C 2,1 (R x [0, +oo)) satisfying 



T 



d 2 F 



d 2 t 



O TT 1 

s as < oo a.s., 



the Ito type formula 

F(B? ,t) = F(0,0)+J*^F(B?,s)dB? 

(2.2) + jf ^F(Bf, s)ds + hJ* ^F(B*, s)s 2H - 1 ds 

holds. For the Ito formula of general Gaussian processes, for example, see Alos et al pQ 
and Nualart-Taqqu [22]. 

Recall that fBm B H has a local time C H (x, t) continuous in (x, t) E R x [0, oo) which 
satisfies the occupation formula (see Geman-Horowitz [9]) 

(2.3) I c/)(Bf,s)ds= [ dx [ c/)(x,s)C H (x,ds) 
Jo Jr Jo 

for every continuous and bounded function (f>(x, t) : R X R + — * R, and such that 



C H (x,t) 



rt i 

/ 8(Bf - x)ds = lim— \(s € [Q,t], \Bf - x\ < e), 
Jo £ J-° 2e 



where A denotes Lebesgue measure and 5(x) is the Dirac delta function. Define the 
so-called weig hted local time S£ H (x, t) of B H at x as follows 

Sf H (x, t) = 2H f s 2H - l C H {x, ds) = 2H f 5{B^ - x)s 2H - 1 ds. 
Jo Jo 

Then the the occupation formula ()2.3f) can be rewritten as 

(2.4) 2H I <f>(B? ,s)s 2H - 1 ds= I dx I <f>(x, s)^ H {x, ds), 

Jo Jr Jo 

and the following Tanaka formula holds: 

(2.5) \B? -x\ = \x\+ [ sign(#f - x)dB? + ££ H {x, t). 

Jo 

Consider the integral representation of fBm B H of the form 

(2.6) Bf = [ K H (t,u)dB u , 0<t<T, 

Jo 
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where B is a standard Brownian motion and the kernel Kff(t,u) satisfies 

( , 7) „)_ KH (I (t . uf -i 



with a normalizing constant kh > given by 

2HTa - H) 



3 ,« X V2 



r(# + i)r(2 - 2fT) y 

Define the operator on L 2 ([0;T]) by 

r H , T h(t) = J' ^K H (t, u)h(u)du, h G L 2 ([0; T]), 

and let S denote the set of all smooth functions on [0, T] with bounded derivatives. 
Then the function Yn,Th{t) is continuous and the transpose T* H t of Tjj,t restricted to 
the interval [0, t] (0 < t < T) is given by 

T* HT g{u) = —khu^~~ h — f m H ~i(m — u) H ~^g(m)dm, 0<u<t 

du J u 

for g G S, the set of all smooth functions on [0, T] with bounded derivatives. In 
particular, for ^ < H < 1 we have 

1 1 3 

T* Ht g(u) = (H )khuz~ h / m H ~2(m — u) H ~2g(m)dm. 

The following lemma gives the computation of Tn.t^Ht ( see 03] > f° r ^ ne P ro °f)- 
Lemma 2.1. Le£ / G S. Then, for ^ < H < 1 we have 

(2.8) r^rj^/Ca) = fT(2ff - 1) / 

J o 

Finally, recall that fractional Malliavin derivative for smooth random variable 

F = f {Io m '"'So r]n{t)dBtH ) ' 

with / G C fe °°(R n ) is defined 

(2.9) DfF = E|£ (j\(t)dB?> • - , f\ n (t)dB^ rub), 

where 77^- = 1,2,..., n) satisfies J Q T T* HT r]j(t) dt < 00. Denote = Th,t^* ht D^ ' ■ 
We have (Proposition 6.24 in |14j ) 

[ g(s)dB?) =[ (E[T* H . T g(s)]) 2 ds+ [ [ E [B? g(r)B? g(s)] ds 
Jo / Jo Jo Jo 

= a H I [ E[g(s)g(r)]\s-r\ 2H ~ 2 dsdr+ [ [ E [Bf 'g(r)B? g(s)] dsdr 
Jo Jo Jo Jo 

with an = H(2H — 1), if the right hand side of this identity is finite. 



ft 

\s-u\ 2H - 2 f(u)du. 
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3. Power variation of fractional local time 

To use Young integration to establish integral with respect to local time of fBm, we 
first investigate p- variation of the mapping x i— ► Ji? H (x, t) for every t > 0. 

Definition 3.1. Let p ^ 1 be a fixed real number. A function f : [a, b] t— > R is o/ 

bounded p-variation if 



sup^ - < 00, 



where the supremum is taken over all partition A n = {a = xq < x\ < ■ ■ ■ < x n = b} of 
[a,b]. 

By the local nondeterminacy of fBm we can prove the following estimate. 
Lemma 3.1. For all s,r € [0, T], s > r and < H < 1 u>e /iai>e 

(3.1) 8 2H r 2H -v?>K(8-r) 2H r 2H , 
where \i = E(Bf B 1 ^) and k > is a constant. 
Lemma 3.2. For t > 0,x G R sei 



Then the estimate 

(3.2) £7[(5f (6) - 5f (a)) 2 ] < C H>t (& - a) 1+Q 

ZioWs /or all^<H<l,0<a< 2H H ^ 1 and a,b £ M, a < b, where Cjj,t > is a constant 
depending only on H,t, so the process {B^(x) : x € M} has a-Holder continuous paths 
with a e (0, for every 0<t<T. 

Proof. For a,b G R, a < b define the function f(x) = l( a& ](x). A straightforward 
calculation shows that 

Dfl [ o )00 )(B/ r ) = *(B/ r )l [0)t] ( S ). 

It follows that 



E (B?{b) - B?(a)) 2 = E (jf* f{Bf)dBf 

= H(2H -\) f f\ s - r\ 2H - 2 El {aM i(B?,B?)dsdr 
Jo Jo 
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for \ < H < 1. Now let us estimate Elum (Bf, B^). We have 



b ,-b 



E \a,b] 2 { B ^ , B^ 



1 



a J a 27T^ 



exp 



r 2H x 2 - 2fixy + s zn y z ) ) dxdy 



2p 2 



1 fW i. 



br^ — fix 



27T / ^ 



< 



< 



< 



1 fW 



2irJ^ 



r H (b - a) 
P 

r (a-l)H 



(b-a) 



dx I 

J a, 


P 

" — fix 




P 


x ' 2 dx 


( 1 


b 

j: 








|1+Q 





e z y2 dy 



271 J °-r H -fix 



e i y dy 



where p = \J s 2H r 2H — p 2 . It follows from Lemma 13,11 that for ^ < H < 1 

_ _ rt 

E{B?{b)-B?{a)) 2 < AH{2H - l)n(b - a) 1+a 



ds I ( s -r) 2H - aH ~ 2 r- H dr 



4K(2H-l) tH(1 __ a) / r 1 
(1 " «) Wo 



'0 Jo 

(1 - x) 2H - aH - 2 x- H dx ) (6 - a) 1+ " 



This obtains the estimate (|3.2p for | < H < 1. 



□ 



Theorem 3.1. Let Jz? H (x,t) be the weighted local time of fBm B H with | < H < 1. 
T/ien /or every < t < T , the limit 

lim Y E\^ H {a i+l ,t)-^ H {a h t)\™-' 



i=0 



exists, and so 



n-l 



lim V |^ H (a m ,t) -£ ,H (a i ,t)\ P = in L 1 

n— >oo ^— ' 1 1 

i=0 

/orp > 3 n^-i j where {ao,ai, ■ ■ ■ ,a n } is a partition of[a,b] such that maxj{|ai+i — Oj|} 
as re tends to infinity. 

Proof. Keeping the notation in Lemma [3.2j. we get 

jSf"(x,t) = 2 - Bf (a) 

by fractional Meyer- Tanaka formula (|2.5p . where 
(3.3) 0t(x) = (Sf -x)+-(-x)+ 

It follows that for any i > 0, x S R, 

Tl - 1 2H 



^|i?"(a m ,t)-if H (a,,t) 



i=0 



<C H 



'n-l 



V |0*(a m ) - &(a0|3»=i + V Bf (ai+i) - Bf (a*) 
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Noting that the function (j>t(x) is Lipschitz continuous in x with Lipschitz constant 2, 
we get the convergence in L 1 

n—1 n—1 

2^ \<pt{a i+ i) - &(ai)|a»-i < C H 2^ \a i+ i - a»| 53=1 < C#| An\ (b - a) — >0, 
i=0 i=0 

as n —* oo. We finally have 

n-l 



SfCoi+O-SfCof) 3 "- 1 <C7 H (6-a) 

by Lemma l3.2i This completes the proof. □ 

Theorem 3.2. Let | < if < 1. Then the weighted local time ^f H (x,t) of fBm B H is 
of bounded p-variation in x for any < t < T , for all p > ^zj, almost surely. 

Proof. Keeping the notation in the proof of Theorem 13.11 we have 

Sf H (t,x)=2(d H (x)-Bf(x)). 

Let [-N, N] contain the support of ^f H (t, x) in x and let 

D := {-N = a < ai < ■ ■ ■ < a n = N} 

be a partition of [— N, N]. Then 

sup V|<M«m) - M^)\ p < sup2f - a t y ^ 2 P (2N) P < oo. 

D i D i 

On the other hand, Lemma 13.21 yields 
B t H (x)-B t H (y) 



l + m 

<G H ,T\x-y\ 2 , a.s. 



for all x, y £ K by Garsia-Rodemich-Rumsey Lemma, where Gh,t is a nonnegative 
random variable such that E(G 9 HT ) < oo for all q > 1, and < m < 2H ^ 1 ■ It follows 
that for all p > ^ > ^ 



sup^(Bf (a i+ i) - Bf(ai)) p < G^ T sup ^(a m - a*) 



p(l + m) 
2 



D ; D 



„ p(l + m) 

= G P HT {2N)^^ < oo, a.s. 
Thus, we complete the proof. □ 

Corollary 3.1. Let ^ < H < 1. T/ien i/te usual local time £ H (x,t) is of bounded 
p-variation in x for any < t < T and all p > , almost surely. 

Now we can establish one parameter integral of local times of fractional Brownian 
motion. Recall that the Riemann-Stieltjes integral (see Young |27j ) 



f{x)dg{x) := lim ^ f(£ i+1 )(g(x i+ i) - gfa)) 
|A„|->0 . =Q 



exists if / and g have finite p-variation and finite q- variation in the interval [a, b] respec- 
tively, and / and g have no common discontinuities, where € [xi,Xi + i], p, q ^ 1, 
| + | > 1, | A„| = maxosjisjn - »t|- 
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Proposition 3.1. For ^ < H < 1, if fix) is of bounded p-variation in x with 1 ^ p < 
Yrg, then the Young integral 

f f(x)JZ> H (dx,t) 

J a 

exists in L 1 for any < t < T . 

Consider the properties and characterization of one parameter integral. We first have 



„ n 

/ f(x)J? H (dx,t)= lim y2f(x j . 1 ){^ H (x j ,t)-^ H (x j ^,t)) 

Jr An ^ J^{ 

n n— 1 

Y,ffa-i)& H fas) -j2ffa)^ H fa^) 



lim 



3=1 
n 



3=0 



= ~ A Um nE Ufa) ~ ffa-l))^ H fa,t) 

by adding some points in the partition A n to make 

Sf H (x u t) = 0, & H (x n ,t) = 0, 
which yields the following 

Corollary 3.2. Under the conditions of Proposition [3A\ the integral 

^ H (x,t)df(x) 

K 

exists in L 1 for any < t < T and 

f(x)Sf H {dx,t) = - [ Sf H (x,t)df(x). 
Jr 

In particular, for f £ C 1 (R) we have 

(3.4) / f(x)^ H (dx,t) = - [ ^ H (x,t)f'(x)dx. 

Jr Jr 

Define the modifier p by 

cexp( (j; _ 1 1 ^_ 1 ), x 6(0,2), 



(3.5) 



p{x) 



0. 



x 0(0,2), 



where c is a normalizing constant such that J^p(x)dx = 1. Set p n (x) = np(nx). For a 
locally integrable function g(x) we define 



9n(x) = I p n (x - y)g(y)dy, n > 1. 
Then g n {x) is smooth and 



g n {x) = / p(z)g(x )dz, n > 1. 

Jo n 

Then, similar to the proof of Theorem 2.1 in [7j one can gives the following 
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Lemma 3.3. Let g(x) be a measurable function with bounded p-variation, where 1 ^ 
p < j^j- Suppose that g n (x) is defined as above, then we have 

lim / g n (x)£? H (dx,t) = / g(x)£f H (dx,t) a.s, 
Jr Jr 



Recall that if F is the difference of two convex functions (This means that F is an 
absolutely continuous function with derivative of bounded variation.), then the Ito- 
Tanaka formula 

F(Sf) = F(0)+ f F'_{B?)dBf + \ ! J? H (x,t)F"(dx) 

JO z Jr 



F(0)+ f F'_{Bf)dB* -~ [ F'_(x)J? H (dx,t) 
Jo 1 Jr 



holds. This is given by Coutin et al [I] (see also Hu et al [13]). The following result is 
a modification of this formula. 

Proposition 3.2. Let the measurable function f : R — * R be of bounded p-variation 
with 1 ^ p < Y—H an< ^ ^ F ^ e an absolutely continuous function with derivative 
F' = f. If the integral f(B^)dB^ exists, then 

F(B?) = F(0)+ f f{B*)dBf- \ I f(x)J? H (dx,t). 
Jo 1 Jr 

Proof. For n ^ 1 we set 

F n (x) : = / p n {x -y)F{y)dy, 
Jr 

where p n is the mollifier defined in (|3.5p . Then, it is well-known that for each x, 

lim FJx) = F(x), lim F'Jx) = fix). 

n— >oo n— >oo 

Moreover, F n E C 2 (R) and the fractional ltd formula yields 

(3.6) F n (B?)=F n (0)+ f F' n {B»)dB^ + H f s^F^B^ds. 

Jo Jo 

If n tends to infinity, then it is easy to see that F n (B^) converges to F{B^) almost 
surely, and F^{B^) converges to f(B^). 

Next we consider the limit of the last term in (13.61) . We have 



2H f s 2H ~ l F' 1 [{B^)ds= f ^ H (x,t)F^(x)dx 
Jo Jr 

F' n {x)^ H {dx,t) 

f(x)^ H (dx,t) 
almost surely, as n tends to infinity. Finally, we see that 

f F' n {Bf)dBf = f F' n (B*) o W^ds - f f(B?)dB?, in (S)' 
Jo Jo Jo 



and also almost surely as n — > oo, by the convergence of the other terms in f|3.6|) . So, 
the result follows. □ 
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4. Weighted quadratic covariation 
In this section, we use one parameter integral 

f(x)£> H (dx,t) 



to study the weighted quadratic covariation \f(B H ),B H ] . Let B be a standard 
Brownian motion. Recall that if / is locally square integrable, then the quadratic 
covariation [f(B),B] of f(B) and B exists in L 1 and 

[f(B),B] t = - [ f(x)J? B (dx,t), 

JR 

where J£ B {x,t) is the local time of B and [f(B),B] is defined by 



n-l 



[f(B),B] t = lim ^2{f(B (k+1)t/n ) - f(B kt/n )}(B (k +l)t/n ~ B kt / n ) 



n— >oo 

fe=0 



with the limit being uniform in probability. For this see Follmer et al [8] and Eisen- 
baum [5]. However, this is not true for fBm B H , i.e. in general 



[f(B H ),B H ] t ^- [ f(x)J? H (d. 

JR 



'x,t), 0<t<T, 



because [f(B H ),B H ] t = for 1/2 < H < 1. But, as we pointed out before from 
Gradinaru-Nourdin [10], Gradinaru et al [11] and Mishura-Valkeila [19] one can drive 
an inspiration to construct the weighted quadratic covariation of f(B H ) and B H . 

Definition 4.1. Define the process [f(B H ), B H ] J W) , < t < T by 

[f{B H \B H ] { p 

n— 1 

: = 2H ^L E ^{/(^D/J " f(B^ n )}(Bf (k+1)/n - Bf k/n ) 



n—roo ■ 

k=0 



with the limit being uniform in probability. This process is called the weighted quadratic 
covariation of f(B H ) and B . 

Clearly, if H = ^ the weighted quadratic covariation coincides with the usual qua- 
dratic covariation of Brownian motion B, and in particular we have 

[B H ,B H ]f ] =t 2H . 

A classical result (see Klein and Gine |15] when p = 1) shows that, for all t £ [0, T] 

n-l 

A n (T, t) := 2H £ k 2H -H m (kT/n)(B? k+1)T/n - B^J* 



k=0 



= 2H Yl ^ 1 ( B kt/(n±) ~ B kt/(n±)) 
k=0 

— > H2 P t 2Hp , a.s. 

as n tends to infinity, and moreover, A n (T,t) also converges to ^2pt 2Hp in L 2 , where 
a.2 P is the 2p- moment of a standard Gaussian random variable G ~ -/V(0, 1). Thus, the 
following theorem is also due to Gradinaru-Nourdin [10J. 
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Theorem 4.1. Let g G C(R) and let p > 1 a integer. Then the convergence 

(4.1) £ k™^g{B» t/n ) [Bf k+l)t/n - — > M2p / giB?)****-^ a.s. 

/jo/ds /or all < H < 1, as n tends to infinity. 

As a direct consequence of the above theorem, for all / G C 1 (1R) we see that 

Urn £fc 2H -H/(5& + i)/n) " f(B? k/n )}(B? (k+1)/n - B* /n ) — ff{Bf 
almost surely, as n — > 00, and 

(4.2) [f(B H ),B H ]^ = - f f(x)J? H (dx,t) 

Jr 

by occupation formula and the identity (|3.4p . More generally, we have 

Theorem 4.2. Let the measurable function x 1— > f(x) be of bounded p-variation with 
1 ^ p < -tth- Then the weighted quadratic covariation \f(B H ),B H ]^ W ^ exists, and we 
have 

(4.3) [f(B H ),B H ]W = - f f(x)J? H (dx,t) 

Jr 

for all t G [0,T]. 

Proof. If / G C 1 (R), this is Let now / £ C 1 (E). For m ^ 1 we set 



f m (x) := / p m {x -y)f(y)dy, 
Jr 

where p m is the mollifier defined as in (13. 5p . Then / m G C 1 (]R) and we have 



[f m {B H ),B H ] ( p = - [ f m (x)J? H (dx,t) 

J R 



for all t G [0, T]. Moreover using Lebesgue's dominated convergence theorem, one can 
prove that as m — > 00, for each x, f m (x) — > f(x). In order to finish this proof we 
consider the double sequence 

ra-l 

a mn (t) := £ k 2H ~ l {fm(Bf k+1)t/n ) - f m {B* /n )) (Bf k+1)t/n - < /n ), m,n > 1. 

fc=0 

Elementary analysis show that the following convergence hold almost surely: 
hm a mn {t) = [f m (B H ),B H } { ™ ] = - [ f m (x)J? H (dx,t), m > 1, 

lim [f m {B H ),B H ] { p =- f f(x)J? H (dx,t), 
n-1 

hm cw(t) = £ A: 2 *" 1 {/(Bg +1)t/n ) - (^ +1)t/n - Bg /B ), n > 1, 

fc=0 

where the first convergence above is uniform in m. These imply that 

lim lim a mn (t) = lim lim a mn (t) 

m— >oo n— >oo n— »oo m-^oo 

for all £ > 0, and theorem follows. □ 
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According to Proposition 13.21 we get an extension of Ito's formula. 

Corollary 4.1. Let the measurable function f : R — > R be of bounded p-variation with 
1 ^ p < y^tj an d let F be an absolutely continuous function with the derivative F' = f . 
If the integral f(B^)dB^ exists, then we have 

(4.4) F(B?) = F(0) + jT f(B*)dB? + \ [f(B H ), B H )f } 

for all t G [0,T]. 



Clearly, this is an analogue of Follmer-Protter-Shiryayev's formula (see [8]). It is 
an improvement in terms of the hypothesis on / and it is also quite interesting itself. 
Recall that a process X is called a finite quadratic variation process if its quadratic 
variation [X, X] exists. FBm with Hurst index H > ^ is of finite quadratic variation. 
If X is a finite quadratic variation process and if F G C 2 (R), then the following Ito's 
formula holds: 

(4.5) F{X t ) = F(0) + J* f(X s )d~X s + i [f(X),X] t , 

where f = F' and the integral Jq f(X s )d~X s is the forward (pathwise) integral defined 
by 

t ^ ft 



f{X s )d-X s = lim - / f(X s )(X s+£ - X 8 )ds 
o £ l° £ Jo 

and [f(X),X] t = f* f'(X s )d[X, X] s . If X is a continuous semimartingale, then the 

integral 

f f(X s )d~X s = [ f{X s )dX s 
Jo Jo 

is classical Ito's integral. For fBm B H with | < H < 1 we have 

ff{Bf)d-B*= ff{Bf)dB? + H f f\B»)s^ds. 
Jo Jo Jo 

However, the formula (|4,5p is only effective on twice-differentiable functions. Our results 

is a modification of this formula, but our Ito formula is only effective on fBm with 

| < H < 1. Thus, it is natural to ask whether the ltd formula (|4.5p for any absolutely 

continuous function 



F(x)=F(0)+ f f(y)dy 
Jo 



holds or not. This is an interesting question. We refer to Russo and Vallois [25j (see 
also Gradinaru et al [HI [12]) for a complete description of stochastic calculus with 
respect to finite quadratic variation process. 

At the end of this section, we consider the process Y H of the form 

Y} H) := f*f(Bf)dBf + [f{B H ),B H ]f ) , < t < T. 
Jo 

Clearly, when H = ^, the process l^ 1 / 2 ) is the backward integral J*f(B s )d*B s with 
respect to Brownian motion B = {Bt, < t < 1}, defined by 

/ f(B s )d*B s := lim £ f(B u+1 )(B u+1 - B u ) 

JU 0<ti<t 2 <--"<tn<t 
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in probability, where A n = maxj{|ij — tj_i|}. Follmer-Protter-Shiryayev have proved 
in [8 J (see also Eisenbaum [5]) that if / is locally square integrable, then the above limit 
in probability exists, and 



[ f(B s )d*B s = - [ f(B s )dB s 
Jo Ji-t 



io Ji-t 
for t € [0, 1], where Bt = -Bi-j is the time reversal on [0, 1] of the Brownian motion B. 

Theorem 4.3. Let the measurable function x t— > f(x) be of bounded p-variation with 
1 < y^j- If the integral Jq f{Bf)dBf exists, then we have 



- [ f(B H s )dB H s = [ f(B?)dB? + [f(B H ),B H ] { p 

JT-t JO 

for all t 6 [0, T], and so, 

[ f(x)J? H (dx,t) = f f(B?)dB?+ [ T f{B S s )dB~S S) 

Jr JO JT-t 

where B H t = B^_ t is the time reversal on [0, T] of fBm B H . 
Proof If / G C 1 (R) , then the Ito formula for F with F' = f yields 

F{B~H t )=F{B~H Q )+ f f{B~H s )dB~H s + H f f '{B* S )(T - sf^ds 
Jo Jo 

for all t £ [0, T]. Elementary calculus yields 

[ T f(B~H s )dBH s = F(0) - F(B?) — H f T f'(B~H s )(T - sf^ds 
JT-t JT-t 

= - f f(B?)dB?-[f(B"),B H ]^. 
Jo 



For more details about the Ito formula of general Gaussian processes, see Alos et al [lj 
and Nualart-Taqqu [22J. 
Now, for / C X (R) we set 



fn{x)= / p n {x - y)f(y)dy, n > 1, 
Jr 

where p n (x) = np{nx) and p(x) is the mollifier defined in (|3.5p . Then f n (x) is smooth 
and 

fn(x) > f(x) 

for all tends to infinity, and 

- f f n (B^ s )dB^ s = f f n (B?)dB?+[f n (B H ),B H }^ ] 

JT-t Jo 

= f f n (Bf)dBf- [ f n (x)J? H (dx,t) 
Jo Jr 

for all n > 1. This completes the proof by taking n — > oo. □ 
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5. TWO PARAMETER INTEGRALS OF LOCAL TIME 

Now we turn to the local time-space integral. We will use some idea from Feng- 
Zhao [7j to define the two parameter integrals of local time 

rt r 

g(x,s)^f H (dx,ds). 



Recall that the function (x,y) i— > F(x,y), defined on [a,b] x [c, d] is of bounded p- 
variation in x uniformly in y, if 

m 

sup \F(xi,y) - F(xi-i,y)\ p < do, 

2/G[c,d] i=l 

where {o = xq < x\ < ■ ■ ■ < x m = b} is an arbitrary partition of [a, b], and furthermore, 
it is of bounded p, q- variation in (x, y), if 

n / m \ 1 

su p ^i Ai? (^'%)i p <°°> 

[a,b]x[c,d] j=1 \ i=1 J 

where 

AF(xi,yj) = F(xi,yj) - F(x i - 1 ,y j ) - F(xi,yj-i) + F(» i _ 1 ,y i _i), 

and {a = xq < x\ < ■ • • < x m = b; c = yo < y% < • • • < y n = d} is an arbitrary partition 
of [a, b] x [c, d]. 

A function (x, y) i— ► /(x, y) is called to have a jump at (xi, yi) if there exists an 5 > 
such that for any e > 0, there exists (x2, 2/2) satisfying max{|xi — x%\, \yi — 2/2 1 j- < $ and 

I/O2, 2/2) - f(xi,y 2 ) - f{x 2 , yi) + f(xi, yi)\ > e. 

A function (x,y) 1— ► f(x,y) is called to satisfy t/ie finite large jump condition if for 
any e > 0, there exists at most finite many points {x±, . . . , x ni }, {yi, ■ ■ ■ ,y mi } and 
a constant 5 £ > such that the total (p, g)-variation of / on [x, x + 5 £ ] x [y', y"] is 
smaller than e when [x, x + 5 e ] n {xi, . . . , x ni } = 0, and the total (p, g)-variation of / 
on [x', x"\ x [y, y + <5 e ] is smaller than e when [y, y + <5 E ] n {yi, . . . , y m i} = 0- 

Let now G:RxR^tbea continuous function of bounded gi-variation in x 
uniformly in y, and be of bounded ^-variation in y uniformly in x\ F : R x R — > R 
be of bounded p, q- variation in (x, y) and satisfy the finite large jump condition, where 
P,q,qi,q2 > 1. Then the Young integral (see Theorem 3.1 in Feng-Zhao [7j) 

pb pd n m 

/ / G(x,s)F(dx,dy) := ]im y2^2G(xi-i,yj-i)AF(xj,yj) 

J a Jc A m , n ->-0 

is well-defined, where A mi „ = maxjj{|(xj, yj) — (xj_i, y,_i)|}, if there exist two mono- 
tone increasing functions p : R — > R+ and <r : R — ► R+ such that 



E. 



n,m 



1 \ / 1 \ 1 

O" Tl— —T-, T-, < OO. 



7T 



1/91 / \ m 1 / 92 7 n 1 /Pm 1 / p9 



More works for two-parameter p, q- variation path integrals can be fund in Feng-Zhao [TJ. 
By taking p(u) = ti a and o~(u) = with a £ (0, 1), one can prove the following. 
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Proposition 5.1. Let the measurable function F : R x [0, t] — > R be of bounded p,q- 
variation in (x, t), and of bounded ^-variation in x uniformly in t, and satisfy the finite 
large jump condition. If these parameters p, q, 7 satisfy the conditions 

2H 

(5.1) 1<7<^ EF> 2Hpq<2Hq + 3H-l, p,q>l, 

1 — tt 

then the Young integral of two parameters 



Sf H {x,s)F(dx,ds) 



is well-defined, almost surely, and 

F(x,s)^ H (dx,ds) = [ [ ^ H (x,s)F(dx,ds) 
Jo Jr 

^ H {x,t)F{dx,t). 

Corollary 5.1. Under the conditions of Proposition HOI we have 

[ [ F(x,s)^ H (dx,ds) = - I dx I ^-F(x, s)^ H \x,ds), 
Jo Jr Jr Jo " x 

provided F{x,t) G C 1 ' 1 ^ x [0,T]). 

Corollary 5.2. Under the conditions of Proposition HOI we set 

F m (x,s):= / / p(r)p(z)F(x - —,s - —)drdz, m ^ 1, 
Jo Jo mm 

where p is the mollifier defined in (I3.5j) . Then the convergence 

[ [ F m {x,s)5? H (dx,ds) — *l [ F{x,s)^ H {dx,ds) 
Jo Jr Jo Jr 

holds almost surely, as m — > 00. 

By using Corollary 15.21 and occupation formula, we immediately get an extension of 
the Ito formula as follows. 

Theorem 5.1. Let F G C 1 ' 1 ^ x R+). Suppose that the function (x,t) 1— > -§^F(x,t) is 
of bounded p, q-variation in (x, t), and of bounded ^-variation in x uniformly in t, where 
these parameters p,q,j satisfy the condition (|5.ip . If the integral f Q J^F(Bf , s)dBf 
exists, then the following ltd formula holds: 

F(Bf,t) = F(0, 0) + J* §-/(B?, s)ds + J* ^F(Bf, s)dB^ 

'Iff ^-F(x,s)^ H (dx : ds). 
z Jo Jr ax 



Proof. Let F m be defined in Corollary EJ for m > 1. Then F m G C 2,1 (R x R+) and 
d_ 

dx ' 



Jr-F m , m > 1 are of bounded p, q- variation in (x,t). By occupation formula it follows 
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that 

"t pa r rt f)2 



2H J ^F m (B^s)s 2H - 1 ds = J^dx J — Fm (x,s)J? H (x,ds) 

^-F m (x,s)^ H (dx,ds) 
ox 

■^-F(x,s)^ H (dx,ds), 
ox 

almost surly, as m — > oo. On the other hand, the fractional Ito formula yields 
F m (B?,t) = F m (0, 0) + J* ^jF m {Bf, s)ds + J* -^-F m (B^ , s)dB^ 

for all m > 1. This completes the proof by taking m — ► oo. □ 

As the end of this paper, for a measurable function / : R x [0, T] — ► R we consider 
the weighted quadratic covariation \f(B H ,-),B H ] of f(B H ,-) and B H , defined by 

[f(B H ,.),B H ]W 

n—1 



with the limit being uniform in probability for t G [0, T], where {t k := t^,k = 
0, 1, 2, . . . , n} is a special partitions of [0, t]. For Brownian motion B we have known 
that the quadratic covariation 



n — 1 



[f(B,-),B] t = lim J2{f(B tk+1 ,t k+1 ) - f(B tk ,t k )} (B tk+1 - B tk ) 



n— >oo 

k=0 



exists as a limit uniformly in probability, and 

[f(B, •), B\ t = f f(B s ,s)d*B s - I f(B s , s)dB s , t G [0, 1], 
Jo Jo 

if f(-,t) is locally square integrable and f(-,t) is continuous in t as a map from [0, 1] to 
Lf (M.). For this see Eisenbaum [5] and Follmer et al [5]. 

Theorem 5.2. Let f G C(lRxR_|_) be of bounded p, q-variation in (x,t), and of bounded 
^-variation in x uniformly in t. If these parameters p,q,j satisfy the condition (|5.1|) . 
then we have 

(5.2) [f(B H ,.),B H ] ( p = - f f f(x,s)J? H (dx,ds) 

Jo Jm. 

for all t G [0,T]. 

Proof. Let / G C 1,1 (1R x [0, T\). Thanks to occupation formula, it is enough to prove 



(5.3) 



[f(B H , •), B»] f > = 2H jT ±f(B*, s)s^ds 
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for all t E [0, T]. Consider the decomposition 



n-l 



E* 3 * -1 {fK +1 ,t k+1 ) - f(B*,t k )} « +1 - bd 



k=0 

n—1 o n—1 

= E ^If^ - O 2 + E - - <) 

fc=0 fc=0 

+ ~ E **)«, - O + E ^"^x^ - <) 

k=Q k=0 

= A 1 +A 2 + A 3 + Ai. 
According to Theorem 14.11 we have 



k=0 J0 



s)s 2H - 1 ds 



almost surely, as n tends to infinity. Moreover, elementary calculus and the Holder 
continuity of fBm can show that A2,A^,A4 — > a.s., as n tends to infinity. Thus, we 
have obtained f)5.3[) . 

Let now / E" C 1 ' 1 ^ x [0,T]). For m > 1 we define f m as follows 



f 2 f 2 r z 

f m (x,s):= / p(r)p(z)f(x ,s )drdz, m ^ 1, 

Jo Jo mm 

where p is the mollifier defined in ([53]) ■ Then / m E C 1 ' 1 (R x R+) and we have 

[/ m ( J B H ,-),^]r ) =- r / / m (x, S )^(<ix,d S ) 

Jo Jr 

for all i E [0,i\. By applying the convergence of the double sequence 
n-l 

a m n(t) := E {/m« +1 ,4+l) " fm(B? k ,t k )} - m,n > 1, 

fe=0 

we can show that •), 5^] exists and equals to 

f(x,s)5? H (dx,ds). 

This completes the proof. □ 
From the proof of Theorem 15.21 we see that 

[/(** •), B»] t (H/) = 2H jf* A /(B f , s)s ^-l ds 

for all / E C 1 ' 1 (1R x M+). According to Theorem 15.11 and Theorem 15.21 we get an 
analogue of Follmer-Protter-Shiryayev's formula (see [8]). 
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Corollary 5.3. Let F G C 1 ' 1 (R x R + ),f(x,t) = ^F(x,t) and let the conditions in 
Theorem I5.il hold. Then we have 

F(B?,t) = F(0,0)+ [ f(B?,s)dB? 
Jo 

+ f t ^F(B^s)ds+ l -[f{B H ,.),B 



HI (W) 



for all t G [0,T]. 



Theorem 5.3. Zei / G C(Kxl + ) 6e of bounded p, q-variation in (x,t), and of bounded 
^-variation in x uniformly in t. Assume that the integral J Q f(Bf \s)dBf exists. If 
these parameters p,q,j satisfy the condition (|5.ip . then we have 

- [ T f(B^ s ,T- S )dBS s = [ t f(B^,s)dBf+[f(B H ,.),B H ]^ ) 

JT-t Jo 

for all t G [0, T], where B H t = B^_ t is the time reversal on [0, T] of the fBm B . 
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